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Abstract 

We define Discrete Quasi-Einstein metrics (DQE-metrics) as the critical points 
of discrete total curvature functional on triangulated 3-manifolds. We study DQE- 
metrics by introducing some combinatorial curvature flows. We prove that these 
flows produce solutions which converge to discrete quasi-Einstein metrics when the 
initial energy is small enough. The proof relies on a careful analysis of discrete dual- 
Laplacians which we interpret as the Jacobian matrix of the curvature map. As a 
consequence, combinatorial curvature flow provides an algorithm to compute discrete 
sphere packing metrics with prescribed curvatures. 

1 Introduction 

In [¥] , Cooper and Rivin introduced a type of combinatorial scalar curvature which is 
defined as angle defect of solid angles. Their curvature is a combinatorial analogy of 
the scalar curvature in the infinitesimal sense and is closely related to Regge's calculus. 
They used total curvature functional to study the combinatorial analogue of conformal 
deformation of sphere packing metrics. They found that the Hessian of total curvature 
functional is of great importance. 

In [6], Professor Glickenstein introduced a combinatorial version of Yamabe flow. 
The flow is a system of equations connects CR-curvature (Cooper-Rivin's curvature) 
to sphere packing metrics. The evolution of CR-curvature was then derived and shown 
to satisfy a combinatorial heat equation. The discrete Laplacian was defined so as to 
describe the combinatorial heat equation. Moreover, [H| studied the dual structure 
in Euclidean triangulations coming from sphere packing, and expressed the discrete 
Laplacian by the geometric elements in the dual structure. 

In this paper, we will define DQE-metric (discrete quasi-Einstein metric) as the 
critical point of total curvature functional and give some analytical conditions for 
the existence of DQE-metrics. We will study DQE-metric by combinatorial curva- 
ture flow methods. We will mainly concern two types of discrete energies and their 



1 



negative gradient flows. One is the total curvature functional, it determines a combi- 
natorial CR-curvature flow of second order. Another is the discrete quadratic energy, 
a 3-dimensional analogue of discrete Calabi energy introduced in [5] , it determines a 
combinatorial CR-curvature flow of fourth order. These flows exhibit good existence 
and convergence properties, that is, the solutions of these flows converge to discrete 
quasi-Einstein metrics when the initial normalized quadratic energy is small enough. 
Thus these flows provide useful tools for the study of manifolds with sphere pack- 
ing metrics. It's interesting that we can design algorithm to compute discrete sphere 
packing metrics with user prescribed curvatures because these flows can be modified 
as user's will so as the modified flows tend to find user prescribed target curvatures au- 
tomatically. It's also interesting that much of our process can be formally generalized 
to higher dimensions and to other definitions of curvatures. 



2 Conformal structure and DQE- metric 

2.1 Sphere packing metric and edge-tangential-sphere 

Consider a compact manifold M of dimension 3 with a triangulation T on M. The 
triangulation is written as T = {V, E 7 F 7 T}, where the symbols V,E 7 F 7 T represent 
the set of vertices, edges, faces and tetrahedrons respectively. Throughout this paper, 
all functions defined on vertices are regarded as a column vector and N = V#, the 
number of vertices. Furthermore, all vertices marked by v\, vn are supposed to be 
ordered one by one. We often write i instead of v\. 

Sphere packing metric is a map r : V — > (0, +oo) such that each edge {i,j} € E 
between vertices i and j has length lij = r< + Vj, and each tetrahedron {i,j, k,l} G T 
with six edge lengths Uj, Uk, hi, Ijk, lju Iki can be realized as an Euclidean or hyperbolic 
tetrahedron. In [4], Cooper- Rivin called tetrahedrons generated by this way "confor- 
mal" . Mathematicians had found that a tetrahedron is conformal iff there exists a 
unique sphere tangent to all of the edges of the tetrahedron. Moreover, the point of 
tangency with the edge {i,j} is at distance r» from Vi. 

We may think of sphere packing metrics as points in II > , N times Cartesian 
product of (0,oo). However, not all points in R> represent sphere packing metrics, 
because we still need some nondegenerate conditions besides positive condition. In 
fact, consider a conformal tetrahedron k,l} £ T generated by four positive radii 
Ti , rj , rfc , ti , it can be realized as an Euclidean tetrahedron of same edge lengths if and 
only if 




so, the space of admissible sphere packing metrics is 

DJl r = { r e R^o I Qijki > 0, V{i,j, k,l}eT}. 

In [3], Cooper and Rivin showed that 9Jl r is an open set in R> ' simplex connected, 
but not convex. 
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Remark 1. For a n-simplex, the conformal condition, i.e. the existence of sphere 
packing metric, is equivalent to the existence of a n dimensional sphere tangent to 
each of the edges. 



2.2 CR-curvature and total curvature functional 

Given a Euclidean tetrahedron {i,j, k, 1} € T, the solid angle of vertex i is denoted by 
ctijki ■ Cooper and Rivin's combinatorial scalar curvature is defined to be angle defect 
of solid angles. Concretely, the curvature Ki at a vertex i is 

Ki = 4-7T - OLijM, 

{i,j,k,l}eT 

here, the sum is taken over all tetrahedrons having i as one of its vertex. We write this 
curvature as "CR-curvature" for short. The curvature Ki is a combinatorial analogue 
of scalar curvature in smooth cases, for it measures the difference between, on the one 
hand, the growth rate of a small ball centered at vertex {i} in M and, on the other 
hand, the growth rate of the volume of an Euclidean ball of the same radius. 

Cooper and Rivin first defined this type of combinatorial curvature in [4], they 
also studied the sphere packing metric structure and proved the metric cannot be 
deformed while keeping the solid angles fixed. Their central conclusion were derived 
from a careful study of the following total curvature functional 

N 
i=l 

They showed that the functional is weakly concave and is strongly concave up to 
scaling. In fact, by the famous Schlafli formula, one can get 

V r S = K, 

d{K x ,--- ,K N ) 



Hess r S = 



d(n,--- ,r N ) 

From [3], [Tl], [7] or Lemma [4.41 we proved in this paper, the curvature map K = 
K(r) = V r S is locally an embedding when restricting to the space transverse to 
the scaling direction r. This means the combinatorial scalar curvature can locally 
determine the sphere packing metric up to a scaling. 

Glickenstein introduced the so called combinatorial Yamabe flow 

dri v 
dt 

to study combinatorial Yamabe problem in jB]. This geometric flow connects CR- 
curvatures and sphere packing metrics. Glickenstein showed that the evolution of 
CR-curvature satisfies a heat equation driven by discrete dual-Laplacian. 

[5], [3] and [B] provided great perspectives and deep inspirations for the study 
of combinatorial manifolds with circle packing metrics. The combinatorial structure 
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of the triangulation consists of finite vertices, edges, faces, tetrahedrons and their 
adjacent relations. We can always express all the finite relations by incidence matrix, 
so methodologically, it's more convenient to use the language of matrix when dealing 
problems arisen in combinatorial manifolds. Under this belief, we will give a new 
interpretation of discrete dual-Laplacian as the Jacobian matrix of the curvature map. 



3 Discrete quasi-Einstein metric 

Finding good metrics is always a central topic in geometry. For a smooth manifold, 
Einstein metric is a good candidate for privileged metric on the manifold. On an 
Einstein manifold, the Riemannian metric g is proportionate to the Ricci curvature, 
that is Ric = Xg. 

For a combinatorial manifold (M,T), consider sphere packing metric r as an ana- 
logue of Riemannian metric g, while CR-curvature K as Ricci curvature Ric, we write 
K = Xr formally as the analogue of Einstein condition Ric = Xg. 

Definition 3.1. Given (M 3 ,T), where M 3 is a compact manifold, T is a fixed tri- 
angulation on M 3 . A sphere packing metric r with K — Xr is called a DQE-metric 
(discrete quasi-Einstein metric). A DQE-metric K = Xr is called "flat", "positive" 
or "negative" when the constant A satisfies A = 0, A>0orA<0 respectively. 

Notice that, when r is a DQE-metric, i.e. K — Xr, then A = jpjp-- Besides 
their formal similarity to smooth Einstein metrics, DQE-metrics themselves have great 
importance. They are the critical points of the normalized total curvature functional. 

Theorem 3.2. Consider the r-normalized total curvature functional S T = j^t- 

(1) If t = 1, the critical points of S T are DQE-metrics, vice versa. 

(2) If r 1, the critical points of S T are fiat DQE-metrics, vice versa. 
Proof. 

VrS T = j^(V lS ^r) = ^{K *Lr). 

r is the critical point of S T if and only if K — j|7p r '- Thus the conclusions in the 
theorem follow from S = r T K = tS. □ 
We want to know whether DQE-metric exists on a fixed (M, T). Inspired by 
Thurston's work in [15], we believe that there are combinatorial and topological 
obstructions for the existence of DQE-metrics. As combinatorial curvature flow 
methods ([2], [6]) provide efficient tools in the study of circle packing metrics, 
we will study DQE-metrics by introducing some combinatorial CR-curvature 
flows which are negative gradient flows of some discrete energies in the following 
sections. 



4 



We mainly concern two types of discrete energies and their negative gradient 
flows. One is the total curvature functional S = r T K, it determines a combi- 
natorial CR-curvature flow of second order. Another is the discrete quadratic 
energy C = ||i^|| 2 , it determines a combinatorial CR-curvature flow of fourth 
order. The latter seems more powerful than the former, so we introduce the 
fourth order flow first. 



4 Combinatorial CR-curvature flow of fourth order 
4.1 Definition 



Consider the discrete quadratic energy functional: 



N 



C(r) = \\Kf = J2K] 



(4.1) 



Choose coordinate transformations Ui = lnrj, it's a homeomorphism from 
space 9Jt r to space We will use r-coordinate or u-coordinate freely and not 
distinguish the quadratic energy functional as a function of r or as a function of 
u, people can distinguish them from the context without confusion. Remember, 
all variables r, u, K defined on vertices are considered as iV-dimensional column 
vector. 

Under u-coordinate, the Jacobian matrix of CR-curvature map K = K{u) 

is 

/ 

dui 

d(K u ---,K N ) 



ijjNxN 



d(ui, • • • ,u N ) 



9K N 



Ml \ 

duff 



8K N 



then we have 



/ 9C_ \ 
/ <9«i \ 



l dC I \ ur\i 

\ dufj / \ Sum 



dKi 


dK N 


dui 


dui 


dKi 


dK N 


du jv 





J 



\ Kn J 



2L T K. 



Definition 4.1. Combinatorial CR-curvature flow of fourth order is defined 
to be the following gradient flow 



u(t) = --V u C = -L T K, 



(4.2) 
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where "•" means time derivative, and the derivative of column vector means 
derivative of its components. 

The CR-curvature evolves according to K = —LL T K. LL T is a discrete 
fourth order differential operator which acts on function (defined on vertices) 
by matrix multiplication, this is why we call flow (|4.2p of fourth order. 

We call L discrete dual-Laplacian and write "DDL" for short. It is some- 
what abrupt to give L this name, however, L is exactly a type of discrete 
Laplace operator whose definition relies on the dual structure of sphere pack- 
ing metric and conformal tetrahedrons. Next we will look at L intensively and 
interpret why the name "DDL" came about. 



4.2 Discrete dual-Laplacian 

Regard the 1-skeleton of the triangulation as a 3 dimensional graph. In spectral 
graph theory, discrete Laplacians are always taking in the following form (see 

131) 

\f. Y, ii.jiij in. 

where Bij are weights along each edge {i, j}. "j ~ i" means two vertices i and j 
are adjacent, and the sum is taken over every edge having i as one of its vertex. 
Weights function B a , gives each edge j ~ i a weight B^, different weights 
function give different types of discrete Laplacians. Conformal tetrahedron 
has a strong dual structure, this dual structure determines a special weights 
function and the corresponding DDL. 

For adjacent vertices i ~ j, the (i, j)-entry of L is 



dKi d w2{i,j,k,l}£T a ijkl 



diij 



E 



da 



ijkl 



dv,j 

{i,j,k,l}eT 



E 



da. 



ijkl 



{i,j,k,i}er 



dr-i 



Let 



Bi 



E 



da 



ijkl 



{i,j,k,i}eT 



dr,\ 



Lemma 4.2. For any 1 < i, j < N, 



L 



i.) 



B ik , 


j = i 


kr^i 






j ~ i 


o, 


else 



(4.3) 
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Figure 1: Dual structure in conformal tetrahedron 

Proof. Solid angles are homogeneous functions of degree zero of sphere pack- 
ing metrics. Use Euler's homogeneous function theorem, we get 

doti jk i daijki , da ijk i , da ijk i . , . . s 

n » + rj J + r k J +n J = 0, (4.4) 
dri OTj or k ori 

which implies La = Ba = ^ Bij. □ 

Remark 2. Glickenstein had already found above formula 14.41 using Schlafli 
formula in [6]. 

Corollary 4.3. L(l, • • • , 1) T = 0. 

Given a conformal 3-simplex r = k, I}, let C denote the center of its 
edge-tangential-sphere, E and F are foots of perpendicular from C to face 
{i,j,k} and face {i,j, 1} respectively. The 2-plane spanned by C,E,F meets 
edge {i,j} perpendicularly at point D. Then E (F) is the center of circle 
inscribed in the Euclidean triangle {i,j,k} In a single conformal 3- 

simplex r, the dual j} T to the edge {i, j} (see [6]) is defined to be the region 
surrounded by four points C, E, D, F. The measurement of *{i,j} T is just the 
directed area of this region, which is denoted by Ayki (See Figure [1]). The area 
is negative when the center of edge-tangential-sphere is at the outside of the 
conformal tetrahedron, and positive when at inside. In the whole triangulation 
(M, T) with fixed sphere packing metric, the dual j} is the sum of all such 
*{i,j} T , where r runs over all conformal tetrahedrons containing {i,j} as one 
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edge. Denote l*j as the measurement of then 



Glickenstein proved ([6]) 



l*j — Aijkl- 
{i,j,k,l}eT 

da-ijki Aijki 



Therefore 



It implies that 



{i,jJj}£T dTj nlij 



I*. 

TiBij = — = TjBji 
Hj 

is symmetric in the two indices. 

Glickenstein defined a discrete Laplace operator "A" in [6] and showed that 

Using Lemma 14.21 we obtain 

A/, = Bijifj ./';) I Lufi + Y, L afi ) = (~ L f)i- 

Regard the discrete Laplace operator A as a matrix (it acts on function / by 
matrix multiplication), the above formula is just the components form of the 
following 

A _ L _ d(K lr -- ,K N ) ^ 
d(ui, ■ ■ ■ ,u N ) ' 

Thus the discrete Laplace operator " — A" is interpreted as the Jacobian of 
CR-curvature map. 

Denote R = diag(ri, ...,rjy), L = RL. A = i?A = — L is called weighted 

K U -,K N ) 
(n>— . r iv) 



DDL. Let A = d 2 Cl '"' ' Kn ? , then L = AR and L = RL = RAR. 



Lemma 4.4. The matrix L = RL is positive semi-definite, rank N—l, and the 
nullspace is the span of the vector (1, • • • , 1) T . Moreover, the matrix equation 
L T x = have a unique nonzero solution x = tr up to scaling, where t can be 
any constant. 
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Proof. It's sufficient to prove that A is positive semi-definite. For any confor- 
mal tetrahedron k,l} G T, it had been shown in the appendix of [7] (or 
see [1] and [2]) that the following 4x4 matrix 



A 



ijkl 



d(ai,aj,a k ,ai 
d(ri,rj,r k ,ri) 



TFr~ 



7FF~ 



dai \ 



9a t i 
^7 / 



is positive semi-definite, has rank 3 and the nullspace is the span of the vector 
(r i ,r j ,r k ,r l ) T . 

We want to extend the matrix A^-jy to a N x N matrix. Remember that all 
vertices are ordered and marked by v%, fjv> we ma Y wen suppose V{, Vj, v k , v\ 
arise at i, j, k, I position respectively in the ordered sequence vi, ...,vj\[. Then 
we get a N x N matrix by putting {A.ijki) s t = ~^rf a ^ ^ ne (s,t)-entry position 
for any s,t £ {i,j, k, I}, and putting at other entries. Without confusion, we 
may still write the extended N x N matrix as A«^, then we have 



A 



d(K u --- ,K N ) 
d(n, • • • ,r N ) 



E 

{i,j,k,l}&T 



A 



ijkl- 



Because each component in the sum is positive semi-definite matrixes, A is 
also positive semi-definite . 

If Ax = 0, then Ayfcjx = 0, i.e. x £ Ker(Aij k i) for any {i, j,k, 1} G T. 
Hence there is a constant t^i s.t. (xi,Xj,x k ,xi) = tijki(ri,rj,r k ,ri). Because 
the manifold M is connected, all must be equal. Thus x = tr, which 
implies Ker(A) = {tr\t S R} and rank(A) = N — 1. 

Since Ar = and A is symmetry, then r T A = r T LRr 1 = 0, hence r T L = 0. 
So Lx = has a unique solution x = tr up to scaling. □ 

Due to Lemma 14.41 we may think L has full rank in the space transverse 
to scaling. Also, matrix A is positive definite in the space transverse to r = 
(ri, • • ■ ,rN) T direction, or equivalently, A is positive definite up to scaling. 

Corollary 4.5. Along the combinatorial CR-curvature flow 14.21 the quadratic 



energy functional is descending. Furthermore, both X^=i M « ana - Wi=i r * are 



riV 



constants. 

Remark 3. By Schlafli formula, the differential 1-form co 



N 



i Kidri 



dS 
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is exact. Moreover, because 9JT r is simply connected, we have 

r N 

S= / J2 Kidri + C > 
Ja i=i 

where a is an arbitrarily selected point form 9JT r . 

Remark 4. For any positive constant C > 0, Jli^i r i = C determines a 
hypersurface in *DT r C R^. At each point r of this hypersurface, the direction 
toward which r scales is transverse to the hypersurface ]X =1 rj = C. Thus L 
has full rank, and A is positive definite along the combinatorial CR-curvature 
ffow[OJ 

4.3 Convergence to DQE-metrics 

We had shown that DQE-metrics are critical points of r-normalized total CR- 
curvature functionals. Notice that, DQE-metrics are critical points of the 
discrete quadratic energy C too. Due to Lemma [4,41 V U C = 2L T K = has only 
DQEM solutions, that is K = Xr. Next we will give an equivalent analytical 
condition for the existence of DQE-metrics by combinatorial CR-curvature flow 
methods. 

Theorem 4.6. If the solution to combinatorial CR-curvature flow 14.21 exists 
for all time and converges to a non-degenerate sphere packing metric r(+oo), 
then DQE-metric exists. Moreover, r(+oo) is indeed one. 
Proof. Along combinatorial CR-curvature flow u = —L T K, we have 

L = —LL T K 

and 

C = -2K T LL T K < 0. 

The existence of C(+oo) follows from the existence of r(+oo). This implies 
C(+oo) = due to the fact that C is descending along flow u = —L T K. Hence 
L T (+oo)K(+co) = 0, using Lemma |4~41 we get K(+oo) = Ar(+oo). □ 
We have shown that the existence of DQE-metrics is necessary for the con- 
vergence of combinatorial CR-curvature flow. Actually, it is almost sufficient. 
In fact, we have 

Theorem 4.7. Suppose there is a DQE-metric r* such that K* = X*r*. Then 
the solution of the CR-curvature flow (|4.2p exists for all time and converges to 
r* if the initial metric r(0) is close to r*, or equivalently, the initial normal- 
ized quadratic energy C(r(0)) = \\K (0) - K*\\ 2 = J2iLi( K M ~ K if is sma11 
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enough. 

Proof. Notice that the combinatorial CR-curvature flow (|4.2p is an au- 
tonomous ODE system u = T(u) = —L T K. We want to prove u* is a local 
attractor of the flow. By calculation we find 

r(u*) = —L T K* = -cL T r = 0, 

and 

D u *T(u) = -L T L < 0. 

Although —L T L is only negative semi-definite in the whole space 9Jt r , it is 
negative definite along the curvature flow due to the transversality announced 
in Remark [H Therefore, u* is a local attractor of the flow. The system is 
asymptotically stable at u*. The following four conditions 

• The initial metric r(0) is close to r* 

• The initial metric u(0) is close to u* 

• The initial CR-curvature K(0) is close to K* 

• The initial normalized quadratic energy metric C(0) is small enough 

are all equivalent due to the locally diffeomorphism between K and u. They 
imply both the long time existence and convergence of combinatorial CR- 
curvature flow (I4.2j) . □ 

Remark 5. Remember A = LR~ l = g lf 11 "' ' Kn ? is positive semi-definite, we 
can also consider combinatorial curvature flow r = — ^V r C = —AK. This 
flow is also a fourth order combinatorial CR-curvature flow, it has completely 
similar properties with combinatorial CR-curvature flow (|4.2p . 



5 Combinatorial CR-curvature flow of second order 

Definition 5.1. The second order combinatorial CR-curvature flow is 

r = -K. (5.1) 

We also consider the normalized flow 

r = \r-K (5.2) 

where A = impr- Since = 2r T r = 2r T (Xr — K) = 0, then ||r|| 2 is a constant 
along the normalized flow (|5.2p . By calculation, we get K = —AK. A is a 
second order discrete differential operator, it's why we call flow (|5.2p of second 
order. Moreover, we have C = -2K T AK, S = -XS - C = -\\K - Xr\\ 2 < 0. 
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Theorem 5.2. If the solution of combinatorial CR-curvature flow (|5.2p exists 
for all time and converges to a non-degenerate sphere packing metric r(+oo), 
then DQE-metric exists. Moreover, r(+oo) is indeed one. 
Proof. Along the normalized flow ([53]) . S = -XS - C = -\\K - Xr\\ 2 < 0. 
Under the hypothesis that r(t) — > r(+oo) we know S(+oo) = 0, which implies 
K(+oo) = Ar(+oo). □ 
Let T(r) = \r — K = — K, we have 

D r T(r) = —D r K + XD r r + r(D r X) T 

. .Ar + K 2S , T 

= -A + XI N + r( - j. jjrr) 

= XI N - A + - 



rr T rQg - Ar) T 
X{1n - ir-^) -AH — — . 



At DQE-metric point r* we have 

T 

L» r r| r »=A*(J JV -^)-A. (5.3) 

Next we give two sufficient conditions so as the normalized flow (|5.2p exists 
for all time and converges to DQE-metric. 

Theorem 5.3. Given (M 3 , T). Let r* be a DQE-metric with K* = X*r*,X* < 
0. Assuming the initial normalized quadratic energy C(r(0)) = ||-f^(0) — K*\\ 2 
is small enough, then the the solution of the normalized flow (|5.2p exists for 
t £ [0, +oo) and converges to r*. 

Proof. Notice that In — rr T /\\r\\ 2 is always positive semi-definite. As A* < 
at r*, then 

T 

D r T\ r , =X*(I N -^)-A<0. 

Moreover, rank (D r T\ rt ) = N — 1 and the kernel of -D r r| r< , is exactly tr, t S R. 
Along the normalized flow (|5.2p . ||r|| is a constant, r never scales, we may think 
LVr]^ as a negative positive definite matrix, hence each flat or negative DQE- 
metric r* is a local attractor of the normalized flow (15.21). □ 



Theorem 5.4. Given (M 3 ,T). Let r* be a DQE-metric with K* = X*r*. As- 
suming Ai(A), the first eigenvalue of A at r*, satisfies Ai(A) > A*. Furthermore, 
assuming the initial normalized quadratic energy C(r(0)) = ||i^(0) — K*\\ 2 is 
small enough, then the the solution of the normalized flow (|5.2p exists for 
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t S [0, +00) and converges to r*. 

Proof. Select a orthogonal matrix P, such that P T AP = diag{0, X%, ■ ■ ■ , Ajv_i} 
Write Q = (eo,ei,-- - , ejv-i), where ej is the (i + l)-column of Q. Then 
Aeo = and Aei = A^, 1 < i < N — 1, which implies eo = and 
r _L 1 < i < N — 1. Then (ijv — |prp)ej = e^, 1 < i < N — 1, which implies 

P T {In ~ wv) p = diag{0, 1, • • • ,1}. Therefore, 

T 

AT| r , = X*(I N - „) - A = P T dm 5 {0, A* - Ai,- • • ,A*-Ajv_!}P. 

If Ai(A) > A*, then D r T\ r , < 0, ran/c {D r T\ r ,) = N - 1. The proof can be 
finished by completely the same reason as explained in the proof of Theorem 
IQl □ 



6 Combinatorial G-Laplacian and G-curvature flow 

Besides CR-curvature Ki, there are other types of combinatorial scalar cur- 
vatures. In [10], Glickenstein considered the curvature Ci = K^i. From now 
on, we call Ci "G-curvature" for convenient. CR-curvature is invariant under 
metric scaling, while G-curvature scales in the appropriate way. G-curvature 
is like scalar curvature times the volume measure, and so scales like length. 

The total curvature functional becomes S = YliLi Ki r i = J2iLi Define 
a new discrete quadratic energy 

N 

£ = \\Cf = ^C?. 
1=1 

Now we consider its negative gradient flow 

«(*) = ~\^u£ = —G T C (6.1) 

Here C = (Ci, • • • , C N ) T , G = V u £ = 2G T C. Therefore 

V»5 = C, 
Hess u S = G. 

The matrix G is symmetric since it is the Hessian of the total scalar cur- 
vature functional (or see pH]). We call it G-Laplacian. Remember that 
R = diag{r lr -- ,r N }, A = ^^'.'.f^ , denote Q = diag{Ci,--- ,C N }, then 
by calculation we obtain 

G = n + Z = n + RAR. 
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Lemma 6.1. If Cj > for any i and there is at least one C\ > 0, then the 
matrix G is positive definite. 

Proof. The matrix fi and RAR are both semi-positive definite. The kernel 
of RAR is constant column vector 1 = (1, • • • , 1) T up to scaling, however, this 
vector makes 1 T 01 > 0. □ 
When G is positive definite, good properties will come out locally. Let 

m* r = { r £ Wl r | G(r) > 0} . 

All metrics in this open set exactly determine positive definite G-Laplacians. 
We use to denote it's corresponding domain under coordinate transforma- 
tions Ui = lnr,. 

For any convex set U C 9It*, since G = Hess u S is positive definite, the 
G-curvature map C(u) = V U S is an embedding on U. This means locally 
every G-curvature determines a unique metric when G-curvature is positive 
definite. G-curvature contains the r term, hence is not invariant under scaling. 
This local rigidity property of G-curvature takes scaling into account, and is 
different with the case of Cooper and Rivin's curvature. 

Now we give an example which is an analogue of shrinking sphere in the 
case of smooth Ricci flow. 

Example 1. Given (M 3 ,T), and r(0) is the initial sphere packing metric. 
Let C(0) be the corresponding initial G-curvature and G(0) the initial G- 
Laplacian. If G(0)C(0) is a constant column vector, then the combinatorial 
G-curvature flow has a shrinking or steady solution. Suppose G(0)C(0) = Al, 
since AiV = 1 T G(0)C*(0) = l T ft(0)C(0) = C(0) T C(0) = £(0), we get 

Then the G-curvature flow ii = —G T C with initial metric u(0) has a unique 
solution 

r(0) 1 

In fact we can solve the above G-curvature flow equation by method of 
undetermined coefficients. The solution is shrinking when A > and steady 
when A = 0. Notice that, the condition that G(0)C(0) is a constant column 
vector can be interpreted as "isotropy", it is an analogue of the standard 
"round" sphere in a certain sense. □ 

It is conceivable the G-curvature flow will exhibit good local convergence 
properties when the G-Laplacian is positive definite. For any user prescribed 
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target G-curvature with positive G-Laplacians, the modified G-curvature flow 
tends to find the target G-curvature automatically. 

Theorem 6.2. Given an admissible G-curvature C with positive G-Laplacian, 
that is, there is a metric u £ 50?* that determines C = C(u). Consider the 
modified quadratic energy 

N 

i=i 

The modified combinatorial G-curvature flow 

u{t) = -\^u£ = G T {C-C) (6.2) 

exists for all t > and converges to u as t — > +oo whenever the initial quadratic 
energy is small enough. 

Proof. Notice that the G-Laplacian G = G(u) is positive definite, then use 
the methods in the proof of Theorem 14.71 we can get the above results. □ 

Remark 6. For any user prescribed target G-curvature C (suppose C is well 
selected, i.e. it is admissible and it determines positive definite G-Laplacian), 
any algorithm aiming to minimize £ = Xyi=i(Ci — Ci) 2 , such as Newton's 
gradient descent algorithm, tends to find sphere packing metric f automatically. 

7 Hyperbolic geometric background 

For hyperbolic manifolds, it's more natural to built the triangulation by hyper- 
bolic tetrahedrons, each block is embedded in H 3 . For a hyperbolic conformal 
tetrahedron {i,j, k,l} £ T, Cooper-Rivin showed that the matrix 

d(ai,aj,a k ,ai) 
a{n,rj,r k ,ri) 

is positive definite in [3] and [II) . So 

{i,j,k,i}eT 

is also positive definite due to the same reason expressed in the proof of Lemma 
14.41 Notice that the total curvature functional is changed to 

N 

S = 2vol{M) + s ^K i r i . 

i=i 
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Using Schlafli formula again, we have 



N / N 



) 



N 



dS = Y^ K i dr i + 2dvol(M) +^2ndKi 




and hence 

V T S = K, Hess r S = A. 

Therefore, the curvature map K = K(r) is also an embedding and hence locally 
a diffeomorphism. This implies the following 

Theorem 7.1. In hyperbolic geometric background, the sphere packing metric 
which determines constant zero curvature is isolated. □ 
Consider the following three discrete curvature flows 

• r = —K. 

• r = -\V r C = -AK, 

• u = -±V„C = —L T K. 

Here L = AR^, = diag {sinh(ri), • • • , sinh(rTv)} and the coordinate change 
is Ui = Intanh^f. It's interesting that these three discrete curvature flows 
possess similar convergence properties, we only list the results here without 
proof. 

Theorem 7.2. If the solution to any of the above three curvature flows exists 
for all time and converges to a non-degenerate metric r(+oo), then there exists 
constant zero metric and r(+oo) is indeed one. 

Theorem 7.3. Suppose there is a constant zero metric r* , then these three 
discrete curvature flows exist for all time and converge to r* when the initial 
energy is sufficiently small. 

Remark 7. As was discussed in section [6l we can also investigate G-curvature 
flow in hyperbolic geometric background. In this case, the G-curvature should 
be Ci = Ki sinhrj. We omit the details since the processes and the conclusions 
are almost the same as before. 



8 Future Work 



We pointed out something more we want to do in the future. 
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8.1 Find combinatorial condition so as DQE-metric exists 

We do not even know whether DQE-metrics exist for a triangulated compact 
manifold M 3 . Inspired by Thurston's work in [15] (or see [2]), there are topo- 
logical and combinatorial obstruction for the existence of constant curvature 
metric, we want to seek similar topological and combinatorial condition for 
the existence of DQE-metrics. We also want to know how many DQE-metrics 
there are in a conformal class. 

8.2 How to define discrete curvatures in high dimensions? Dose 
discrete curvatures carry any topological information? 

Consider a n-dimensional triangulated manifold (M, T) with a given flat cone 
metric, such that it becomes a PL-manifold. If n = 2, the only singularity 
arises at vertices. However, in the case n > 3, singularities arise both at 
vertices, edges and even (n — 2)-dimensional faces. What is the "right" analog 
of the sectional curvature, Ricci curvature and scalar curvature? It seems the 
"right" definition of curvatures is variant according to different purpose, such 
as classification purpose, capturing topological information purpose, finding 
good metric purpose and so on, and there are not solely one "right" definition 
of curvature. 

8.3 How to estimate Ai(A)? 

We proved that the DQE-metrics r* with Ai(A) > A* is a local attractor of 
the second order combinatorial CR-curvature flow (15. 2D . We want to develop 
efficient methods to estimate the lower bound of Ai(A). We also want to know 
if there are topological and combinatorial obstructions for Ai(A) > A*. 

8.4 How far can the "flow" method go? 

We still consider a ra-dimensional triangulated manifold (M, T) ■ Suppose we 
have defined a type of scalar curvature K : V — > H for each sphere packing 
metric r. Then K determines a discrete energy C(r) = Ei=i K l Similarly, we 
consider the negative gradient flow of this energy 

u = -^V U C = —L T K. 

Then 

C = —2K T LL T K < 0. 
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Notice that L = ... and rank(LL ) = rank(L). If the combinatorial 
scalar curvature K defined to be have the property rank(L) = N — 1, then 
this flow converges to Ker(L T ) when the initial energy is sufficiently close to 
Ker(L T ). 

The properties of the matrix L is of most importance. The matrix L came 
from the relation between curvature K and the metric r (or u), The full rank 
(up to scaling) of L is equivalent to locally rigidity, i.e. curvature determines 
metric. If a curvature K is well defined and can be used for the purpose of 
classifying PL-manifolds, then it should determine the PL-metric in a certain 
sense. As explained above, "classification", "rigidity", "full rank of L up to 
scaling" are closely related. 

8.5 Combinatorial Gauss-Bonnet formula in higher dimensions 

Suppose the discrete scalar curvatures defined in a certain manner have the 
properties that they are homogeneous functions of degree zero. Then they are 
invariant under scaling of metrics. Three are N vertices and N curvatures 
K\, Kn, however, the essentially independent number of metrics ri,...,r/v 
is N — 1 when modular scaling. Thus the curvatures K\, Kn must have 
some relations, K lies in a no more than (N — l)-dimensional submanifold in 
H N . If the curvatures possess locally rigidity additionally, they must lie in a 
hypersurface in H N . For the case dim(M) = 2, people have already known the 
combinatorial Gauss-Bonnet formula 

N 

Y J Ki = 2v X {M). 
i=i 

Is there an analogue of combinatorial Gauss-Bonnet formula in higher dimen- 
sions? Can we write the hypersurface out clearly? 

8.6 How to get global rigidity? 

The space of admissible sphere packing metrics 9Jt r is not convex, which is the 
main obstruction to get global rigidity. In [T5] . Feng Luo suggested finding 
a good parametrization so as the admissible space of metrics became convex 
in the new coordinate. In Cooper-Rivin considered a new coordinate sys- 
tem ki = The admissible space of metrics became convex under this 
coordinate change. Unfortunately, negativeness of the Hessian matrix of total 
curvature functional under new coordinate seems hard to preserve. It seems 
not easy to seek for a better parametrization so that the admissible space 
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of metrics becomes convex and the positiveness (negativeness) of the Hessian 
matrix is preserved simultaneously. We still need to think up more efficient 
methods to get global properties. 
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